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Abstract 

We give a simple proof to derive the optimal flux which minimizes 
the ground state energy in one dimensional Hubbard model, provided the 
number of particles is even. 

KEY WORDS: Hubbard model, flux phase problem. 

1 Introduction 

We consider the Hubbard model on the ring (i.e., one dimensional system 
with periodic boundary condition), where the magnetic flux is threaded through 
the ring. Our problem is to obtain the optimal flux which minimizes the ground 
state energy. 

To be precise, we define the Hubbard Hamiltonian as follows: 

L L 

H ■= ^ *x,s+i4+i,<r c a l g + (h.c.) + U x n x ,fn x ,i, 
°"=T,1 x=1 x=1 

where L (L > 3) is the number of sites, the site L + 1 is equivalent to the site 
1, t x>x+1 E C, |tx,a;+i| 0, U x e R, 4,a ( res P- °x,<t) is the creation (resp. 
annihilation) operator which satisfies the canonical anticommutation relations, 
and Tlx, (7 • — u^x,it' 

We write t x , x+1 = |ta;, x +i| exp[^ X)X+1 ], 9 XiX+ i £ [0,2tt). Then, the flux 
which penetrates the ring is defined to be <p := 53x=i The ground state 

energy E (in some fixed number of particles N e ) can be regarded as a function 
of ip (and hence we write E = E(ip)), because it does not depend on any choice 
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of {$x,x+i}x=i which satisfies X^ x =i @x,x+i = ¥■ Our aim is to obtain the flux 
f = Popt which attains min^o^) E(tp). We call ip opt as the optimal flux. In 
general, (p opt is not unique, and we will not discuss the uniqueness question in 
this paper. 

There are some closely related problems in the literature (our problem is the 
same as mentioned in (3) below). (1) it appears in a theory of superconductiv- 
ity [AM, W], (2) In the study of the persistent current [FMSWH, K, YF, FK], 
they discussed whether the response of the Hubbard ring to the external field is 
diamagnetic or paramagnetic, and the influence of the electron-electron interac- 
tion to this property, (3) In high dimensional lattice, the flux phase conjecture 
[HLRW] says that the optimal flux per plaquette is equal to the particle density 
per site. This implies that the diamagnetic feature, which widely holds in the 
one particle system, becomes opposite in high electron density regime. This 
conjecture was rigorously proved by Lieb [L] at half filling. Macris-Nachtergaele 
[MN] gave an improved proof of [L]. 

As for the rigorous study of the Hubbard ring (of even length), Lieb- Loss 
[LL] considered free electron case (U x = 0) at half filling, and computed ip opt 
in general situation so that translation invariance is not assumed. They also 
considered what have more complicated geometry such as tree of ring, ladder, 
etc. Lieb-Nachtergaele [LN] computed ip opt also at half filling when U x = U is 
any constant. In this paper, we obtain <p opt when U x and L are arbitrary, while 
N e is even. Due to the hole-particle symmetry, it suffices to consider N e < L. 

Theorem 

Let N e (< L) be even. 

(1) Assume U x < +oo for allx. E(<p) is minimized if ip = (N e /2 + l)ir (mod 
2tt) (resp. (p = N e ir/2) when L is even (resp. L is odd). 

(2) When U x = oo for all x, E(<p) is minimized if ip — 0, w. 

Remarks. 

(I) We can derive the optimal flux in S z ^ subspaces. 

(a) U x < +oo: the optimal flux takes and ir alternatively as S z varies. For 
instance, when N e = 4n, and L is even, then ip op t — ir (S z — 0,2,4, • • •), and 
Part =0 (S z = 1,3,5, •••)■ 

(b) U x = oo: let m :— N^/Ni (N^ (resp. N±) is the number of up (resp. 
down) spins). We suppose iVj > N± here. When m ^ N, <p opt — 2kn/N e , 
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k G Z (in this case, particles can also be regarded as hard core bosons). When 
m e N, </? opt = 2/c7t/(to + 1) — — l)7r (if (m + 1)L is even), and </5 opt = 
{2k - l)n/(m + 1) - (iV e - 1)tt (if (m + l)L is odd), fc e Z. 

(2) When £/ x = +oo, the proof of Theorem tells us that the ground state 
energy is periodic w.r.t. ip with period n (when S z — 0), and period 2n/N e 
(when m £ Z). This fact and its implications are discussed by Kusmartsev and 
Yu-Fowlcr[K, YF]. 

(3) When N e is odd and U x = +oo, we can still derive the optimal flux, and 
the result is the same as stated in Remark (1). 

On the other hand, when U x < +oo, and N e = L (half- filling), we believe 
<fiopt = 7r/2, 37t/2 as some examples imply {e.g., take t x ,x+i- constant and 
U x = 0). However, in general cases, p> pt could be different depending on the 
value of U x . For example, let L = 4, N e = 3, and t x , x +i = t. When U x = 0, 
E{(p) is minimized if and only if ip = ±4 arcsin(l/ y/E) , while in case of U x = +oo, 
E{(p) is minimized if and only if ip = 0, 2-7r/3, 47r/3. 

(4) SU{2) invariance as well as translation invariance is not necessary to 
prove Theorem. We can let t XtX +i — t x x+1 {a =|, J.) depend also on spin 
variable. In this case, our theorem mentions the optimal flux in the S z = 
subspace only. Besides, our Hamiltonian can include the one body potential 
term as well. 

(5) The argument in the proof, together with that in [LM] yields the ground 
state is unique and has spin zero: S — 0, provided the flux p takes the value 
as stated in Theorem. Moreover, if we let E{S) denote the ground state energy 
in spin S subspace, then we have E{S) < E{S + 2). It becomes equality when 
U x = co, ip = 7r (resp. ip = 0), and L: even (resp. L: odd). 

(6) In general, the spin of the ground state is sensitive to the flux. For 
example, let L be even, N e = An + 2, t x ^ x+ \ = t, and U x = co. Then, one can 
show that, (a) when <p = 0, there is a singlet ground state, but no ferromagnetic 
ones, (b) when p = n, there is a ferromagnetic ground state. 

(7) When U x = co, not for all x, the argument of the proof says the following: 
if (t {a; : U x = oo} < L — N e /2, then the result is the same as in Theorem (1). 
Otherwise, the result is the same as in Theorem (2). 

(8) When U x = 0, and i X)X +i = t, E{<p) is maximized if and only if ip = 
N e n/2 (resp. <p = {N e /2 + l)ir), if L is even (resp. L is odd), which should be 
compared with the fact that E{0) — E{ir) when U x = co (Theorem (2)). 



3 



(9) When we let L large, \E(0) - E(ip)\ will behave as 0(1/ L) [LN]. 



In section 2, we give the proof of Theorem, which is very simple. Our problem 
is reduced to consider an one-particle Hamiltonian 7i((p) on the graph G which 
is composed of the basis of iV e -fermion Hilbert space. Theorem follows from 
obtaining the optimal flux of H(<p) on G, by the usual diamagnetic inequality 
argument. 

2 Proof of Theorem 

At first, we consider the case in which L is even and N e = An. Due to the 
5*17(2) invariance, it is sufficient to work on S z = subspace (i.e., JVf = N± = 
2n). We fix the basis of the Hilbert space of iV e -fermions: 



that is, to arrange particles in increasing order w.r.t. the space coordinates. 
Our problem is equivalent to consider the one-particle Hamiltonian 7i(ip) on 
the graph G whose sites are composed of B. 



where s xy (ip) :=< x\H\y >, x, y € B. Two sites x, y £ G are connected by a 
bond if and only if s xy (ip) ^ (we note |s XJ( (<p)| does not depend on <p). For 
given <p e [0, 2n), we fix some {6 x , x +i}%=i such that Y^x=i ®x,x+i = <p, and thus 
we suppose s xy (ip), and hence TL(ip), is determined by ip. Then, it is not hard 
to show that: (1) every circuit in G has even length (because L is even), (2) the 
fluxes in these circuits are always integer multiple of ip '■— <P + ^ n7r + (4n — l)ir. 
In fact, let C be the set of circuits in G which have minimal length. Every 
elements of C is given by fixing all particles which have down (resp. up) spins 
and moving each up (resp. down) spins all together until each spins come to 
next spin. To make it clear, we write down an element of C when L = N e = 4: 



B : = i 4i,<7i4 a ,o i , • • • C x Ne ,a Ne l vac > : 

x\ < x 2 < ■ ■ ■ < x Na , Oi =T,| ; i = 1,- • ■ ,-We}, 



(H(ip)u)(x) := ^2s xy (<p)u(y), 
yets 




vac > 



vac > 



T 

C 2,T C 2,i C 4,T C 4,l l VaC ^ 
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The second term 2nir in the definition of ip comes from the fact that up spins 
jump down spins 2n times on the above process, and each jump causes to put 
(—1) on the corresponding s xy (<p). The third term (An — l)n in the definition 
of ip comes from the fact that the 2n-th up spin jumps all the other (An — 1) 
particles when it moves from the site L to the site 1, because we set the basis 
such that particles are arranged in increasing order. 

On the other hand, because of the inequality: J2 X yeB s xy (p)u(x)u(y) > 
~ J2x y eB \ s xy(p)\\ u ( x )\\ u (y)\i we know that the ground state energy is min- 
imized when all off-diagonal elements s xy (p), x ^ y, are non-positive. Let 

(H-u)(x) := -J^yeB^xyiP^iy)- When ip = ( mod 27r )> H(<f) is unitarily 
equivalent to H-, because the fluxes of all circuits in G are all the same [LL, 
Lemma 2.1]. -0 = (mod 2ir) yields p = ir (mod 2ir). This concludes the proof 
when L is even and N e — An. 

When L is even and N e = An + 2, the only thing we have to do is to replace 
ip in the above argument by ip' := p + (2n + l)ir + (An + l)n. When L is odd, 
then ip (or ip' in case of N e = An + 2), should satisfy ip = tt (mod 2n) to have 
optimal flux, because the minimal length of the circuits in G is odd, and so, the 
flux of 7i_ on every elements of C is 7r. When U x = oo, the minimal length of 
circuits in G is 2L, whose flux is 2<p + 2(N e — l)ir = 2p (mod 2n). □ 

Remarks 

(1) As an alternative proof, one can compute the partition function P(p) :— 
Tr[exp(— p3H)] by using the path integral representation [AL], and show that 
P(ip) is maximized if p takes the value stated in Theorem. This approach has 
been done by [GMMU], where they derived the optimal flux in the Falicov- 
Kimball model. 

(2) When the number of electrons is odd, the fluxes of elements of C are 
different from each other, depending on which spins move in the circuit. For 
example, let L = N e = 2n + 1, = n, and N± = n + 1. By the hole-particle 
transformation only for down spins, we can suppose — N± — n, but now the 
flux of down spins is tt — ip (this situation is similar to that discussed in [FK]). 

Our supposition is the following: the "contribution" to the ground state 
energy from C would cancel each other, and an important contribution would 
come from those circuits where up spins and down spins move together in the 
opposite direction which has flux p — (tt — cp) = 2p — w, and has length 2n in 
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G (the meaning of "contribution" could be clear if we consider Tr[exp(— f3H)] 
instead of the ground state energy). 2ip — -k = would give the minimizing 
energy. However, this supposition would not be easy to prove. 

(3) The proof above relies on the special nature of the ring geometry: there 
is always fixed number of particles on only one loop, so that all circuits on the 
graph G favor the same flux or 7r, depending on cases. However, on more 
complicated systems such as two dimensional lattice, the graph G has so many 
different circuits which favor different fluxes so that our argument does not work 
even if U x = oo, except the Nagaoka-case (N e = |A| — 1, U x = oo [N,T]), where 
the optimal flux is zero everywhere. 

3 Conclusion 

In this paper, we derived the optimal flux <p opt in the Hubbard model on the 
ring. Our result is true in general situation so that the translation invariance is 
not necessary to assume, except the number of particles must be even. In this 
section, we briefly discuss the physical interpretation of our result. 

The result (1) of our theorem is consistent with that of [FMSWH], where 
it is shown that, at half-filling, the current response of the ground state is 
paramagnetic (resp. diamagnetic) when N e = An (resp. An + 2) by numerical 
computation. However, these are not equivalent, especially when iV e = An. In 
fact, [FMSWH] showed, when L — 6, N e = A, and U x > 0, the ground state 
is diamagnetic (this also implies why it is not easy to seek (p which maximizes 
E(ip)). Therefore, our contribution may be that there would be no effects of 
spatial disorder. 

The result (2) of our theorem and Remark (2) after that is already found and 
discussed by [K, YF]. However, our proof gives a different picture: the graph G 
consists of rings of larger lengths, for U = oo prohibits the exchange of particles. 

Finally, our argument gives a ring version of the Lieb-Mattis theorem [LM] 
when ip — ip opt (Remark (5) after Theorem). 
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